Dual-role transcription factors stabilize intermediate expression levels.

Precise control of gene expression levels is essential for normal cell functions, yet how they are defined and tightly maintained, particularly at intermediate levels, remains elusive. Here, using a series of newly developed sequencing, imaging, and functional assays, we uncover a class of transcription factors with dual roles as activators and repressors, referred to as condensate-forming level-regulating dual-action transcription factors (TFs). They reduce high expression but increase low expression to achieve stable intermediate levels. Dual-action TFs directly exert activating and repressing functions via condensate-forming domains that compartmentalize core transcriptional unit selectively. Clinically relevant mutations in these domains, which are linked to a range of developmental disorders, impair condensate selectivity and dual-action TF activity. These results collectively address a fundamental question in expression regulation and demonstrate the potential of level-regulating dual-action TFs as powerful effectors for engineering controlled expression levels.

Adeno-to-squamous transition drives resistance to KRAS inhibition in LKB1 mutant lung cancer.

KRASG12C inhibitors (adagrasib and sotorasib) have shown clinical promise in targeting KRASG12C-mutated lung cancers; however, most patients eventually develop resistance. In lung patients with adenocarcinoma with KRASG12C and STK11/LKB1 co-mutations, we find an enrichment of the squamous cell carcinoma gene signature in pre-treatment biopsies correlates with a poor response to adagrasib. Studies of Lkb1-deficient KRASG12C and KrasG12D lung cancer mouse models and organoids treated with KRAS inhibitors reveal tumors invoke a lineage plasticity program, adeno-to-squamous transition (AST), that enables resistance to KRAS inhibition. Transcriptomic and epigenomic analyses reveal ΔNp63 drives AST and modulates response to KRAS inhibition. We identify an intermediate high-plastic cell state marked by expression of an AST plasticity signature and Krt6a. Notably, expression of the AST plasticity signature and KRT6A at baseline correlates with poor adagrasib responses. These data indicate the role of AST in KRAS inhibitor resistance and provide predictive biomarkers for KRAS-targeted therapies in lung cancer.

ChIPr: accurate prediction of cohesin-mediated 3D genome organization from 2D chromatin features.

The three-dimensional genome organization influences diverse nuclear processes. Here we present Chromatin Interaction Predictor (ChIPr), a suite of regression models based on deep neural networks, random forest, and gradient boosting to predict cohesin-mediated chromatin interaction strength between any two loci in the genome. The predictions of ChIPr correlate well with ChIA-PET data in four cell lines. The standard ChIPr model requires three experimental inputs: ChIP-Seq signals for RAD21, H3K27ac, and H3K27me3 but works well with just RAD21 signal. Integrative analysis reveals novel insights into the role of CTCF motif, its orientation, and CTCF binding on cohesin-mediated chromatin interactions.

Aging-disturbed FUS phase transition impairs hematopoietic stem cells by altering chromatin structure.

ABSTRACT: Aged hematopoietic stem cells (HSCs) exhibit compromised reconstitution capacity. The molecular mechanisms behind this phenomenon are not fully understood. Here, we observed that the expression of FUS is increased in aged HSCs, and enforced FUS recapitulates the phenotype of aged HSCs through arginine-glycine-glycine-mediated aberrant FUS phase transition. By using Fus-gfp mice, we observed that FUShigh HSCs exhibit compromised FUS mobility and resemble aged HSCs both functionally and transcriptionally. The percentage of FUShigh HSCs is increased upon physiological aging and replication stress, and FUSlow HSCs of aged mice exhibit youthful function. Mechanistically, FUShigh HSCs exhibit a different global chromatin organization compared with FUSlow HSCs, which is observed in aged HSCs. Many topologically associating domains (TADs) are merged in aged HSCs because of the compromised binding of CCCTC-binding factor with chromatin, which is invoked by aberrant FUS condensates. It is notable that the transcriptional alteration between FUShigh and FUSlow HSCs originates from the merged TADs and is enriched in HSC aging-related genes. Collectively, this study reveals for the first time that aberrant FUS mobility promotes HSC aging by altering chromatin structure.

SPIRAL: integrating and aligning spatially resolved transcriptomics data across different experiments, conditions, and technologies.

Properly integrating spatially resolved transcriptomics (SRT) generated from different batches into a unified gene-spatial coordinate system could enable the construction of a comprehensive spatial transcriptome atlas. Here, we propose SPIRAL, consisting of two consecutive modules: SPIRAL-integration, with graph domain adaptation-based data integration, and SPIRAL-alignment, with cluster-aware optimal transport-based coordination alignment. We verify SPIRAL with both synthetic and real SRT datasets. By encoding spatial correlations to gene expressions, SPIRAL-integration surpasses state-of-the-art methods in both batch effect removal and joint spatial domain identification. By aligning spots cluster-wise, SPIRAL-alignment achieves more accurate coordinate alignments than existing methods.

Systematic comparison of genome information processing and boundary recognition tools used for genomic island detection.

Genomic islands are fragments of foreign DNA that are found in bacterial and archaeal genomes, and are typically associated with symbiosis or pathogenesis. While numerous genomic island detection methods have been proposed, there has been limited evaluation of the efficiency of the genome information processing and boundary recognition tools. In this study, we conducted a review of the statistical methods involved in genomic signatures, host signature extraction, informative signature selection, divergence measures, and boundary detection steps in genomic island prediction. We compared the performances of these methods on simulated experiments using alien fragments obtained from both artificial and real genomes. Our results indicate that among the nine genomic signatures evaluated, genomic signature frequency and full probability performed the best. However, their performance declined when normalized to their expectations and variances, such as Z-score and composition vector. Based on our experiments of the E. coli genome, we found that the confidence intervals of the window variances achieved the best performance in the signature extraction of the host, with the best confidence interval being 1.5-2 times the standard error. Ordered kurtosis was most effective in selecting informative signatures from a single genome, without requiring prior knowledge from other datasets. Among the three divergence measures evaluated, the two-sample t-test was the most successful, and a non-overlapping window with a small eye window (size 2) was best suited for identifying compositionally distinct regions. Finally, the maximum of the Markovian Jensen-Shannon divergence score, in terms of GC-content bias, was found to make boundary detection faster while maintaining a similar error rate.

Chromatin-associated RNA Dictates the ecDNA Interactome in the Nucleus.

Extrachromosomal DNA (ecDNA) promotes cancer by driving copy number heterogeneity and amplifying oncogenes along with functional enhancers. More recent studies suggest two additional mechanisms for further enhancing their oncogenic potential, one via forming ecDNA hubs to augment oncogene expression 1 and the other through acting as portable enhancers to trans-activate target genes 2. However, it has remained entirely elusive about how ecDNA explores the three-dimensional space of the nucleus and whether different ecDNA have distinct interacting mechanisms. Here, by profiling the DNA-DNA and DNA-RNA interactomes in tumor cells harboring different types of ecDNAs in comparison with similarly amplified homogenously staining regions (HSRs) in the chromosome, we show that specific ecDNA interactome is dictated by ecDNA-borne nascent RNA. We demonstrate that the ecDNA co-amplifying PVT1 and MYC utilize nascent noncoding PVT1 transcripts to mediate specific trans-activation of both ecDNA and chromosomal genes. In contrast, the ecDNA amplifying EGFR is weak in this property because of more efficient splicing to remove chromatin-associated nascent RNA. These findings reveal a noncoding RNA-orchestrated program hijacked by cancer cells to enhance the functional impact of amplified oncogenes and associated regulatory elements.

Enterococcus faecalis Strains with Compromised CRISPR-Cas Defense Emerge under Antibiotic Selection for a CRISPR-Targeted Plasmid.

Enterococcus faecalis is a Gram-positive bacterium that natively colonizes the human gastrointestinal tract and opportunistically causes life-threatening infections. Multidrug-resistant (MDR) E. faecalis strains have emerged that are replete with mobile genetic elements (MGEs). Non-MDR E. faecalis strains frequently possess CRISPR-Cas systems, which reduce the frequency of MGE acquisition. We demonstrated in previous studies that E. faecalis populations can transiently maintain both a functional CRISPR-Cas system and a CRISPR-Cas target. In this study, we used serial passage and deep sequencing to analyze these populations. In the presence of antibiotic selection for the plasmid, mutants with compromised CRISPR-Cas defense and enhanced ability to acquire a second antibiotic resistance plasmid emerged. Conversely, in the absence of selection, the plasmid was lost from wild-type E. faecalis populations but not E. faecalis populations that lacked the cas9 gene. Our results indicate that E. faecalis CRISPR-Cas can become compromised under antibiotic selection, generating populations with enhanced abilities to undergo horizontal gene transfer. IMPORTANCE Enterococcus faecalis is a leading cause of hospital-acquired infections and disseminator of antibiotic resistance plasmids among Gram-positive bacteria. We have previously shown that E. faecalis strains with an active CRISPR-Cas system can prevent plasmid acquisition and thus limit the transmission of antibiotic resistance determinants. However, CRISPR-Cas is not a perfect barrier. In this study, we observed populations of E. faecalis with transient coexistence of CRISPR-Cas and one of its plasmid targets. Our experimental data demonstrate that antibiotic selection results in compromised E. faecalis CRISPR-Cas function, thereby facilitating the acquisition of additional resistance plasmids by E. faecalis.

Aging-induced pseudouridine synthase 10 impairs hematopoietic stem cells.

Aged hematopoietic stem cells (HSC) exhibit compromised reconstitution capacity and differentiation-bias towards myeloid lineage, however, the molecular mechanism behind it remains not fully understood. In this study, we observed that the expression of pseudouridine (Ψ) synthase 10 is increased in aged hematopoietic stem and progenitor cells (HSPC) and enforced protein of Ψ synthase 10 (PUS10) recapitulates the phenotype of aged HSC, which is not achieved by its Ψ synthase activity. Consistently, we observed no difference of transcribed RNA pseudouridylation profile between young and aged HSPC. No significant alteration of hematopoietic homeostasis and HSC function is observed in young Pus10-/- mice, while aged Pus10-/- mice exhibit mild alteration of hematopoietic homeostasis and HSC function. Moreover, we observed that PUS10 is ubiquitinated by E3 ubiquitin ligase CRL4DCAF1 complex and the increase of PUS10 in aged HSPC is due to aging-declined CRL4DCAF1- mediated ubiquitination degradation signaling. Taken together, this study for the first time evaluated the role of PUS10 in HSC aging and function, and provided a novel insight into HSC rejuvenation and its clinical application.

Development and Experimental Validation of a Novel Prognostic Signature for Gastric Cancer.

BACKGROUND: Gastric cancer is a malignant tumor with high morbidity and mortality. Therefore, the accurate recognition of prognostic molecular markers is the key to improving treatment efficacy and prognosis. METHODS: In this study, we developed a stable and robust signature through a series of processes using machine-learning approaches. This PRGS was further experimentally validated in clinical samples and a gastric cancer cell line. RESULTS: The PRGS is an independent risk factor for overall survival that performs reliably and has a robust utility. Notably, PRGS proteins promote cancer cell proliferation by regulating the cell cycle. Besides, the high-risk group displayed a lower tumor purity, higher immune cell infiltration, and lower oncogenic mutation than the low-PRGS group. CONCLUSIONS: This PRGS could be a powerful and robust tool to improve clinical outcomes for individual gastric cancer patients.

SODB facilitates comprehensive exploration of spatial omics data.

Spatial omics technologies generate wealthy but highly complex datasets. Here we present Spatial Omics DataBase (SODB), a web-based platform providing both rich data resources and a suite of interactive data analytical modules. SODB currently maintains >2,400 experiments from >25 spatial omics technologies, which are freely accessible as a unified data format compatible with various computational packages. SODB also provides multiple interactive data analytical modules, especially a unique module, Spatial Omics View (SOView). We conduct comprehensive statistical analyses and illustrate the utility of both basic and advanced analytical modules using multiple spatial omics datasets. We demonstrate SOView utility with brain spatial transcriptomics data and recover known anatomical structures. We further delineate functional tissue domains with associated marker genes that were obscured when analyzed using previous methods. We finally show how SODB may efficiently facilitate computational method development. The SODB website is https://gene.ai.tencent.com/SpatialOmics/ . The command-line package is available at https://pysodb.readthedocs.io/en/latest/ .

SOTIP is a versatile method for microenvironment modeling with spatial omics data.

The rapidly developing spatial omics generated datasets with diverse scales and modalities. However, most existing methods focus on modeling dynamics of single cells while ignore microenvironments (MEs). Here we present SOTIP (Spatial Omics mulTIPle-task analysis), a versatile method incorporating MEs and their interrelationships into a unified graph. Based on this graph, spatial heterogeneity quantification, spatial domain identification, differential microenvironment analysis, and other downstream tasks can be performed. We validate each module's accuracy, robustness, scalability and interpretability on various spatial omics datasets. In two independent mouse cerebral cortex spatial transcriptomics datasets, we reveal a gradient spatial heterogeneity pattern strongly correlated with the cortical depth. In human triple-negative breast cancer spatial proteomics datasets, we identify molecular polarizations and MEs associated with different patient survivals. Overall, by modeling biologically explainable MEs, SOTIP outperforms state-of-art methods and provides some perspectives for spatial omics data exploration and interpretation.

Profiling and characterization of constitutive chromatin-enriched RNAs.

RNA species act as architectural scaffolds for nuclear structures including chromatin in eukaryotic cells. However, the composition and dynamics of tightly bound chromatin-associated RNAs during mitosis remains elusive. Here we report the identification of chromatin-enriched RNA (cheRNAs) by biochemical nuclear fractionation coupled with RNA sequencing in both interphase and mitotic phase of A549 and HeLa-S3 cell lines. We show that highly abundant cheRNAs, mostly small noncoding RNAs, are largely maintained in mitotic chromatin, and constitute a substantial part of chromatin RNA throughout cell cycle. We also show that the mitotic retained cheRNAs tend to be cell type nonspecific and might be involved in chromatin accessibility and epigenetic memory of gene expression control. Therefore, we reveal an unexpected set of cell type-nonspecific mitotic retained chromatin-enriched RNAs. We anticipate that the landscape of RNA composition of chromatin both in interphase and mitotic phase would help understanding structure and function of chromatin.

DAZL regulates proliferation of human primordial germ cells by direct binding to precursor miRNAs and enhances DICER processing activity.

Understanding the molecular and cellular mechanisms of human primordial germ cells (hPGCs) is essential in studying infertility and germ cell tumorigenesis. Many RNA-binding proteins (RBPs) and non-coding RNAs are specifically expressed and functional during hPGC developments. However, the roles and regulatory mechanisms of these RBPs and non-coding RNAs, such as microRNAs (miRNAs), in hPGCs remain elusive. In this study, we reported a new regulatory function of DAZL, a germ cell-specific RBP, in miRNA biogenesis and cell proliferation. First, DAZL co-localized with miRNA let-7a in human PGCs and up-regulated the levels of >100 mature miRNAs, including eight out of nine let-7 family, miR21, miR22, miR125, miR10 and miR199. Purified DAZL directly bound to the loops of precursor miRNAs with sequence specificity of GUU. The binding of DAZL to the precursor miRNA increased the maturation of miRNA by enhancing the cleavage activity of DICER. Furthermore, cell proliferation assay and cell cycle analysis confirmed that DAZL inhibited the proliferation of in vitro PGCs by promoting the maturation of these miRNAs. Evidently, the mature miRNAs up-regulated by DAZL silenced cell proliferation regulators including TRIM71. Moreover, DAZL inhibited germline tumor cell proliferation and teratoma formation. These results demonstrate that DAZL regulates hPGC proliferation by enhancing miRNA processing.

Prostate Cancer Transcriptomic Regulation by the Interplay of Germline Risk Alleles, Somatic Mutations, and 3D Genomic Architecture.

Prostate cancer is one of the most heritable human cancers. Genome-wide association studies have identified at least 185 prostate cancer germline risk alleles, most noncoding. We used integrative three-dimensional (3D) spatial genomics to identify the chromatin interaction targets of 45 prostate cancer risk alleles, 31 of which were associated with the transcriptional regulation of target genes in 565 localized prostate tumors. To supplement these 31, we verified transcriptional targets for 56 additional risk alleles using linear proximity and linkage disequilibrium analysis in localized prostate tumors. Some individual risk alleles influenced multiple target genes; others specifically influenced only distal genes while leaving proximal ones unaffected. Several risk alleles exhibited widespread germline-somatic interactions in transcriptional regulation, having different effects in tumors with loss of PTEN or RB1 relative to those without. These data clarify functional prostate cancer risk alleles in large linkage blocks and outline a strategy to model multidimensional transcriptional regulation. SIGNIFICANCE: Many prostate cancer germline risk alleles are enriched in the noncoding regions of the genome and are hypothesized to regulate transcription. We present a 3D genomics framework to unravel risk SNP function and describe the widespread germline-somatic interplay in transcription control. This article is highlighted in the In This Issue feature, p. 2711.

CTCF functions as an insulator for somatic genes and a chromatin remodeler for pluripotency genes during reprogramming.

CTCF mediates chromatin insulation and long-distance enhancer-promoter (EP) interactions; however, little is known about how these regulatory functions are partitioned among target genes in key biological processes. Here, we show that Ctcf expression is progressively increased during induced pluripotency. In this process, CTCF first functions as a chromatin insulator responsible for direct silencing of the somatic gene expression program and, interestingly, elevated Ctcf expression next ensures chromatin accessibility and contributes to increased EP interactions for a fraction of pluripotency-associated genes. Therefore, CTCF functions in a context-specific manner to modulate the 3D genome to enable cellular reprogramming. We further discover that these context-specific CTCF functions also enlist SMARCA5, an imitation switch (ISWI) chromatin remodeler, together rewiring the epigenome to facilitate cell-fate switch. These findings reveal the dual functions of CTCF in conjunction with a key chromatin remodeler to drive reprogramming toward pluripotency.

Avian influenza viruses suppress innate immunity by inducing trans-transcriptional readthrough via SSU72.

Innate immunity plays critical antiviral roles. The highly virulent avian influenza viruses (AIVs) H5N1, H7N9, and H5N6 can better escape host innate immune responses than the less virulent seasonal H1N1 virus. Here, we report a mechanism by which transcriptional readthrough (TRT)-mediated suppression of innate immunity occurs post AIV infection. By using cell lines, mouse lungs, and patient PBMCs, we showed that genes on the complementary strand ("trans" genes) influenced by TRT were involved in the disruption of host antiviral responses during AIV infection. The trans-TRT enhanced viral lethality, and TRT abolishment increased cell viability and STAT1/2 expression. The viral NS1 protein directly bound to SSU72, and degradation of SSU72 induced TRT. SSU72 overexpression reduced TRT and alleviated mouse lung injury. Our results suggest that AIVs infection induce TRT by reducing SSU72 expression, thereby impairing host immune responses, a molecular mechanism acting through the NS1-SSU72-trans-TRT-STAT1/2 axis. Thus, restoration of SSU72 expression might be a potential strategy for preventing AIV pandemics.

HPVMD-C: a disease-based mutation database of human papillomavirus in China.

Human papillomavirus (HPV) can cause condyloma acuminatum and cervical cancer. Some mutations of these viruses are closely related to the persistent infection of cervical cancer and are ideal cancer vaccine targets. Several databases have been developed to collect HPV sequences, but no HPV mutation database has been published. This paper reports a Chinese HPV mutation database (HPVMD-C), which contains 149 HPV genotypes, 468 HPV mutations, 3409 protein sequences, 4727 domains and 236 epitopes. We analyzed the mutation distribution among HPV genotypes, domains and epitopes. We designed a visualization tool to display these mutations, domains and epitopes and provided more detailed information about the disease, region and related literature. We also proposed an HPV genotype prediction tool, which can predict HPV carcinogenic or non-carcinogenic risk genotypes. We expect that HPVMD-C will complement the existing database and provide valuable resources for HPV vaccine research and cervical cancer treatment. HPVMD-C is freely available at Database URL: http://bioinfo.zstu.edu.cn/hpv.

MyoD is a 3D genome structure organizer for muscle cell identity.

The genome exists as an organized, three-dimensional (3D) dynamic architecture, and each cell type has a unique 3D genome organization that determines its cell identity. An unresolved question is how cell type-specific 3D genome structures are established during development. Here, we analyzed 3D genome structures in muscle cells from mice lacking the muscle lineage transcription factor (TF), MyoD, versus wild-type mice. We show that MyoD functions as a "genome organizer" that specifies 3D genome architecture unique to muscle cell development, and that H3K27ac is insufficient for the establishment of MyoD-induced chromatin loops in muscle cells. Moreover, we present evidence that other cell lineage-specific TFs might also exert functional roles in orchestrating lineage-specific 3D genome organization during development.